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ABSTRACT 

The existence of bright quasars at high redshifts implies that supermassive black holes were 
able to form in the early Universe. Though a number of mechanisms to achieve this have been 
proposed, none yet stands out. A recent suggestion is the formation of quasi-stars, initially 
stellar-mass black holes accreting from hydrostatic giant-like envelopes of gas, formed from 
the monolithic collapse of pre-galactic gas clouds. In this work, we modify the Cambridge 
STARS stellar evolution package to construct detailed models of the evolution of these ob- 
jects. We find that, in all of our models, the black hole inside the envelope is able to reach 
slightly more than one-tenth of the total mass of the system before hydrostatic equilibrium 
breaks down. This breakdown occurs after a few million years of evolution. We show that the 
mechanism which causes the hydrostatic evolution to end is present in polytropic models. We 
also show that the solutions are highly sensitive to the size of the inner boundary radius and 
that no physical solutions exist if the inner boundary is chosen to be less than about 0.3 of the 
Bondi radius. 

Key words: accretion, accretion discs - black hole physics - quasars: general 



1 INTRODUCTION 



Over the last decade, high-redshift surveys have revealed the ex- 
istence of bright quasars at redshift s z > 6 dFan et all 1 20061 ; 
Ijiang et alj 120081 ; IWillott et ail |2010|) . Such observations imply 
the existence of supermassive black holes (SMBHs) with masses 
-Wbh > 10 9 Mq less than 10 9 yr after the Big Bang. It remains 
unclear how these objects became so massive so quickly. Substan- 
tial and ongoing research has led to various possibilities, each sub- 
ject to a number of uncertainties (see the recent review bv lVolonteril 
l20ich . 

One popular possibility is that the SMBHs that power high- 
redshift quasars grew from smaller seed black holes (BHs) that 
were the remnants of the first generation of stars. The Big Bang did 
not produce significant amounts of elements heavier than helium 
so these first stars are expected to have been composed purely of 
hydrogen and helium. They have a number of features which dis- 
tinguish them as a population and are thu s termed population 111 
(pop-Ill) stars (see lBrornrn & Larsonll20Q4l for a review). 

Pop-Ill stars are thought to have formed at the centres of 
smaller dark matter (DM) haloes, with virial temperatures in the 
range 10 3 < T v i r / K < 10 4 , in which the gas cooled through 
molecular hydrogen emission. These stars had masses in the range 
100 < Mj Mq < 1000. Models predict that those with masses 
of more than about 260 Mq undergo pair-instability supernovae 
and leave BH remnants of abo ut half of their progenitor masses 
dFrver. Wooslev & Hegerl 1200 ll) . As the DM haloes continue to 
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merge and grow, so the seed BH accretes gas from its surround- 
ings, settles to the centre of the halo and grows in mass. 

How long would it take such a seed BH to grow to 10 9 Mq? 
If the BH were able to accrete constantly at its Eddington-limited 
rate, then a 100 Mq BH would take about 7 x 10 8 yr to grow into 
a 10 9 Mq SMBH for a typical radiative efficiency of the accre- 
tion (e = 0.1). This is barely the age of the Universe at z = 6 
and already brings the scenario into doubt. In addition, the as- 
sumpt ions of sufficient fuel and constant Eddington accretion are 
weak. Ijohnson & Brorrrrnl |2007) find that the supernova marking 
the formation of a BH causes a d elay of order 100 Myr before the 
BH be gins to accrete efficiently. Milosavli evic. Couch & Bromml 
(2009) find that the accretion on to a BH from a dense protogalac- 
tic cloud is self-limiting and the maximu m accretion rate i s abou t 
32 per cent of the Eddington-limited rate. King & Pringle (2006) 
have shown that, if accretion occurs in small, frequent, randomly- 
aligned episodes, the time constraint is relaxed because the BH spin 
stays low on average, reducing the radiative efficiency. In the case 
of monolithic collapse from a protogalactic cloud, it is not clear that 
the material reaches the centre with random alignments, although 
this may well be the case for subsequent merger-driven accretion. 

It thus appears that pop-Ill seeds cannot grow sufficiently 
rapidly to become the BHs that power high-redshift quasars. How- 
ever, if the first generation of stars pollute the interstellar medium 
with a sufficient metal content, then the second generation of stars 
resembles modern populations. During hierarchical mergers, gas 
builds up in the cores of h aloes and undergoes fra gmentation, 
forming dense stellar clusters dClark. Glover & Klessen 200^). The 
dense environment of these stars can lead to frequent stellar colli- 
sions and thence either directly to a massive BH or to a massive star 
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which leaves a massive BH as its remnant dDevecchi & Volonteril 

12009b . 

ISpolvar. Freese & Gondolol (l2008) suggested that DM annihi- 
lation inside pop-Ill stars can significantly heat them. This could 
delay the star' s arrival on the mai n sequence and allow it to grow to 
a larger mass dFreese et al.ll2008l) . Once the DM is exhausted, nor- 
mal evolution would proceed but larger seed BH masses are possi- 
ble than for typical pop-Ill stars. 

A separate branch of possible paths to early SMBH forma- 
tion stems from the direct collapse of supermassive clouds (M > 
10 5 M ) at the centres of DM haloes with T vir > 10 4 K. If H 2 for- 
mation is suppressed, the only coolant for primordial gas is atomic 
hydrogen. The gas is only able to cool efficiently to about 10 4 K 
dTegmark et al.l 19971) and cosmological simulations show that it 
condenses int o thick, pressure-supported discs that are gravita - 
tionally stable dRegan & Haehneltl2009l ; lwise. Turk & Abell2008l) . 
The cores of such discs are then able to collapse into isolated struc- 
tures with masses that exceed 10 4 Mq. 

The formation of these massive objects hinges on the sup- 
pression of H2 formation, at least until particle densities around 
10 5 cm -3 are reached, whereafter collisional dissociation is 
suffic ient to prevent fragmentation (Schleicher, Spa ans & Gloved 
|2010|) . A popular mechanism for this s uppression is photodisso- 
ciation by an ionizing UV background. Shang, Bryan & Haiman 
(2010) estimate that the necessary specific intensity e xceeds what is 
expect ed on average in the relevant epoch. However, Di ikstra et al.l 
(2008) suggest that the inhomogeneous distribution of ioniz- 
ing sources suppres ses H2 formation in a fraction of haloes. 
Spaa ns & Silkl I2OO6) argue that the self trapping of Lya radiation 
during the collapse keeps the temperature above 10 4 K. The impor- 
tance and effect of H2 remains uncertain but the work cited above 
indicates a number of paths which overcome this hurdle. 

The further central collapse of th e discs described above re - 
quires transport of angular momentum. lLodato & Nataraianl(l2006l) 
claim that, for haloes with T vu /T giis < 1.8, angular momentum 
is transported by local gravitational i nstabilities such that the disc 
mainta ins a state of marginal stability. Begelman, Volo nteri & Reesl 
(2006) instead argue that, when a critical threshold of rotational 
support is achieved, the clouds become unstable to the formation 
of bars, which transport angular momentum outwards and ma- 
terial inwards on dynamical timescales. If the gas remains able 
to cool, then the process can be repeated and a cascade of bars 
is formed. This is the bars-w i thin-b ars mechanism described by 
Shlosman, Frank & Begelman d 19891) and seen in the simulations 
d2008l) . lBegelman & ShlosmarJ d2009h go further to 
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show how it also suppresses fragmentation by maintaining super- 
sonic turbulence and thus relaxes the requirement of the absence of 
H 2 . 

Ultimately, if a large body of gas is able to condense in the 
centre of the DM halo, it can take one of two forms. On the one 
hand it may form a very massive star. If the star is sufficiently mas- 
sive, it is prone to general-relativistic instabili ties and collaps es into 
a BH with 90 per cent of the progenitor mass (Sha pirol2004l) . If the 
total mass of the star exceeds about 10 5 Mq the collapse occurs 
before the end of the main sequence. If the total mass is smaller but 
still exceeds 3.4 x 10 4 M^ , then it collapses during helium burn- 



ing dAppenzeller & Frickefl97ltlBond. Arnett & Carrll984l If the 
total mass is smaller still then the star becomes pair-unstable after 
core helium burning and l eaves a BH of about half the progenitor 
mass dOhkubo et alj2006t) . 

If the rate of mass infall is much higher then the envelope of 
the star does not reach thermal equilibrium during the lifetime of 



the star jBegelmanll2O10h . In this case, the core collapses after hy- 
drogen burning is complete. The structure that remains after core 
collapse is a stellar- mass BH embedded w ithin a giant-like enve- 
lope, or a quasi-star ^Begelman et al. 2006). The attractive feature 
of a quasi-star is that the accretion rate on to the BH is limited 
by the Eddington rate of the entire object, which is initially much 
larger than that of the BH alone. The excess energy is carried away 
by convection. 

Quasi-stars are th e sub ject of this paper. 
iBegelman. Rossi & Armitagd |2008, hereinafter BRA08) ex- 
amined their structure using analytic estimates and basic numerical 
results. We use a fully fledged stellar evolution code to study quasi- 
stars' structure and evolution without many of the assumptions 
made by BRA08. 

In Section[2] we describe the Cambridge STARS code and the 
modifications that we have made. In Section[3] we present the re- 
sults of a fiducial model and make a short theoretical diversion 
to explain some of the model's behaviour in Section|4] We then 
present results with varied parameters in Section|5]and compare our 
models with some of the estimates made by BRA08 in Section [531 
We explore the sensitivity of our models to the inner boundary ra- 
dius in Section|6]before concluding in Section|7] 



2 METHOD 

The models that we report on in this paper were computed with a 
modified versi on of the Cambridge STARS co de. The original code 
was written by Eggle tonl dl97lUl972L 1 19731) and it has been sub - 
stantially updated bv lPols et alj < fl995n and lEldridge & Toutld2004h . 
The principal modification to the code is the change to the interior 
boundary conditions. Normally, stellar evolution codes solve for 
the interior boundary conditions r, m, L r — 0, where r is the ra- 
dial co-ordinate, m is the mass within a radius r and L r is the lumi- 
nosity through the sphere of radius r. We replace these conditions 
with a prescription for the BH's interaction with the surrounding 
gas as described below. 

To derive suitable conditions we presume that the pressure in 
the envelope i s domin ated by the radiation from the accreting BH. 
lLoeb & Rasiol (l994) showed that a radiation-dominated fluid in 
hydrostatic equilibrium, which is not generating energy, must be- 
come convective. Thus the envelope of our quasi-star is approxi- 
mated by a gas with polytropic index n — 3. We derive the bound- 
ary conditions below under this presumption but in the calculations 
we use the adiabatic index (7 = (dlogp/dlog p)s) determined 
self-consistently by the equation-of-state module in the model. 



2.1 Radial boundary condition 

The radius of the inner boundary of the envelope must be the point 
at which some presumption of the code breaks down. In particular, 
we choose a radius within which the gas is no longer expected to be 
in hydrostatic equilibrium. A reasonable choice is the Bondi radius 
tb, the radius at which the thermal energy of the fluid particles 
equals their gravitational potential energy. 



1 2 GmM BH 
— mc, = , 

2 s r B 



so 



2GM BH 



(1) 



(2) 
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where m is the mass of a test particle, c s — \fyp~fp~ is the adiabatic 
sound speed, Af B H is the mass of the BH and G is the gravitational 
constant. Inside this boundary we expect the BH's gravity to over- 
come the thermal motion of the fluid. The average radial velocity 
of the fluid is presumed to be inward, although convective motions 
may result in local motions both towards and away from the BH. 



2.2 Mass boundary condition 

For the mass boundary condition, consider the mass of the gas in- 
side the cavity defined by the Bondi radius, equation ((2). By defi- 
nition 



4:Tvr p(r)dr, 



(3) 



where rs is the Schwarzschild radius. Using a general rela- 
tivistic form of the e quation introduces terms of order rs/rs 
dThorne & Zvtkowlll977t) which we ignore because all our models 
have r s < r B . 

To determine the mass of gas inside the cavity, we must make 
some assumption about the density profile of the material therein 
because the code does not model this region. The quasi-star is sup- 
ported by radiation pressure and is expected to radiate near its Ed- 
dington limit. The Eddington limit for the entire quasi-star is much 
greater than the same limit for the BH alone and excess flux drives 
bulk convective motions. The radial density profile of the accretion 
flow then depends on whether angular momentum is transported 
outward or inward. In the former case, the radial density profile is 

_ 3 

proportional to r 2 whereas, in the latter case, it is proportional 
tor~3 dNaravan et al] : 2000). We presume that the viscosity due to 
small scale magnetic fields is sufficiently large to transport angular 

momentum outwards even if convection transports it inwards and 

_ 3 

thus take pyr) oc r 2 . 

Given the density p{r%) = pa at the inner boundary, the den- 
sity profile must be 



p(r) = Po i^-j ■ (4) 
We evaluate equation (0, presuming rs < rj, to find 

M rav = yporl. (5) 

In a radiation-dominated n = 3 polytrope, the pressure and density 
are related by 



P = 



k 



3(1-/3) 



p 3 = Kp 3 



(6) 



/ \ a/3 4 

(Eddington 1918), where k is Boltzmann's constant, p is the mean 
molecular weight of the gas, ran is the mass of a hydrogen atom 
and /3 = P g /P is the ratio of gas pressure to total pressure. Taking 
the adiabatic sound speed to be c s = \J 4P/3p, evaluating the 
Bondi radius using equation 10 and substituting into equation 10, 
we obtain 



A 



8tt / 3GM B h V 
3 \ 2K ) 



(7) 



which we can use to consider the importance of the mass inside the 
cavity. iFowleil (fi"964) gives P = 4.3(M,/M )-2/ A t, where M, 
is the total mass of the object. For a totally ionized mixture of 70 
per cent hydrogen and 30 per cent helium, p — 0.615. For a total 
quasi-star mass of 10 4 Mq, as in our fiducial results, we find that 



M cav = -Mbh when A4bh ~ 390 Mq. So we must include M cav in 
the mass boundary condition. 

Thus, we set the boundary condition for the mass co-ordinate 

to 



M = Mbh + M ca 



(8) 



where A/ cav is given by equation l[5j. In general, different accretion 
modes give different prescriptions for the density inside the cavity 
and thus different cavity masses. In our fiducial results, we have 
assumed that the density inside the cavity is as described above. 
In Section \5?2\ we construct a model presuming that p(r) a r~5 
and find that changing the density profile inside the cavity has little 
effect. 



2.3 Luminosity boundary condition 

The final boundary condition that must be included is the luminos- 
ity. The luminosity is determined by the mass accretion rate through 
the relationship 



eMc 



(9) 



where c is the speed of light, M is the rate of mass flow across 
the base of the envelope and e is the radiative efficiency, the frac- 
tion of accreted rest mass that is released as energy. This fraction is 
lost from the system as radiation so the total mass of the quasi-star 
decreases over time. The actual rate of accretion on to the BH is 
Mbh = (1 — e)M, the amount of infalling matter less the radi- 
ated energy. The luminosity condition is related to the actual BH 
accretion by 



-Mbhc 2 = e'Mbhc 2 



(10) 



We thus implicitly assume that the material travels from the base of 
the envelope to the event horizon within one timestep. The material 
actually falls inward on a dynamical timescale so this condition is 
already implied by the presumption of hydrostatic equilibrium. 

To specify the mass accretion rate from t he inner boun dary we 
begin with the adiabatic Bondi accretion rate dBondill 19521) . 



M ( GM bh) 2 

AiBon = 47TA C po, 



(ID 



where A c is a factor that d epends on th e adiabatic index 7 as de- 
scribed by equation (18) of lBondilfl952l) . For the case of 7 = 4/3, 
A c = 1/ \/2. Almost all of this flux is carried away from the BH by 
convection. BRA08 point out that the maximum convective flux in 
the material is pc s so that the maximum luminosity is 



Anr B pc s 

4 2 3 
-irr B c s p 

7 

— MsonC 2 . 

7A C 



(12) 
(13) 

(14) 



In order to limit the luminosity to the convective maximum, the ac- 
cretion rate is reduced by a factor c 2 /j\ c e'c 2 . We presume that the 
actual convective flux is some fraction of the maximum computed 
above and thus implement the mass accretion rate 



, n (GM BH ) 2 

e'7 c s c 2 



(15) 



where 77 represents the convective efficiency. For the fiducial run, 
we take 77 = e = 0.1. 
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2.4 Further assumptions 

We conclude this section with a brief discussion of the remaining 
presumptions in the code. The temperatures and densities at the 
base of the envelope are generally too low for meaningful energy 
generation from nuclear reactions to take place, so we did not solve 
the chemical evolution equations in the results presented here. To 
confirm that our approximation is sound, preliminary runs were 
performed with nuclear reactions included. No discernible differ- 
ence was found in the results. 

This does not immediately exclude the possibility of nuclear 
reactions occurring inside the hydrodynamic regio n. Using a tem- 
perature profile T oc r _1 dNaravan et al.l 120000 . we estimated 
the composition changes due to the pp-chains presuming complete 
mixing down to 10 -4 re. We found no significant change to the H 
and He abundances and conclude that the associated energy gener- 
ation is also negligible. Although the temperatures in these regions 
are well over 10 s K, the densities in the region inside rn are typi- 
cally only a few g cm -3 . 

We have not performed detailed calculations that consider the 
CNO-cycle. Prior to its collapse to a BH, the stellar core is expected 
to synthesize sufficient CNO to maint ain hydrostatic equilibrium 
through the hydrogen-burning phase (Be gelmanl 120100 . Most of 
these metals fall on to the black hole unless convection is estab- 
lished with the envelope quickly at core collapse. If the core abun- 
dance is mixed, it is largely diluted by the pristine material in the 
envelope. 

Our model of the interior also neglects any loss of heat via 
neutrino emission. To establish whether this is reasonable, we es- 
timated the total neutrino loss rate using the analytic estimates of 
lltoh et al.1 < fl9960 and integrated them over the interior region for 
the models in the fiducial run. If the flow extends all the way to 
the innermost stable circular orbit (ISCO) of a non-rotating BH, 
nsco = 6rs, we find the neutrino losses are at most 6 per cent of 
the total luminosity. If the flow is truncated at 2nsco, the neutrino 
losses come to less than 0.02 per cent. Such losses would effec- 
tively decrease the radiative efficiency but the accretion rate is prin- 
cipally determined by the convective efficiency. We thus believe our 
model envelopes remain stable against catastrophic neutrino losses. 

If the advection of material across the event horizon is faster 
than the time taken by radiation to diffuse out, then energy can be 
lost to the BH. Such radiation trapping reduces t he support for the 
infalling material and enhances the accretion rate. Beg elmanl dl978l) 
determined an approximate expressio n for this in c rease which is 
supported by the numerical results of iFlammangl J 19841) . We re- 
computed the fiducial run with the appropriate factor included and 
found no difference in the results. 

For now, we have not included non-spherical effects of rota- 
tion inside the hydrodynamic region. Such effects may be impor- 
tant. Convection-dominated accretion flows are expected to have 
moderate angular dependence but do not support mass ejection. 
If convection maintains constant specific angular momentum then 
even Keplerian rotation at a few hundred rs becomes dynamically 
insignificant at rs. We expect any outflows to be impeded by the 
material at or near the base of the envelope and that convective tur- 
bulence preserves approximate spherical symmetry in the vicinity 
of the inner radius. We intend to introduce an approximate treat- 
ment of rotation to the envelope models in future work. 

The STARS code includes a detail ed equation-of-state package 
( Eggl etonetalll973l ; IPols etal.ll 19950 which computes the ioniza- 
tion states of H and He and the contribution of H2. The gas is ev- 
erywhere approximated well by a sum of radiation pressure and the 




Figure 1. Plot of density against radius for models in the fiducial ran with 
M BH / M = 5, 100, 200, 500, 1000 and 1194 (see Table[T}. At the base 
of the envelope the density profile steepens because of the steeper pressure 
gradient required to balance the BH gravity. In the outer layers the density 
is inverted, as discussed by BRA08. Note that the initial model has inner 
radius of 1.66 Rq which is too small to be seen. 



ideal gas law but the ionization state of the material has a substan- 
tial effect on the structure of the envelope. 



3 FIDUCIAL RESULTS 

We begin the exposition of our results by selecting a run which we 
shall use to demonstrate the qualitative features of a quasi-star en- 
velope's structure and evolution. We subsequently experiment with 
various properties of the model to indicate how the behaviour is 
affected by such changes. 

The results presented in this section describe a model quasi- 
star with initial total mass (BH, cavity gas and envelope) M 4 = 
1O 4 M . The BH initially has mass 0.0005M. = 5M but the 
evolution does not depend on this fraction. The gas is uniformly 
composed of 0.7 hydrogen and 0.3 helium by mass. The envelope 
is allowed to relax to thermal equilibrium before the BH begins 
accreting. 



3.1 Structure 

In the model, the luminosity is approximately equal to the Ed- 
dington luminosity at the boundary of the innermost convective 
layer. The accretion rate varies between about 1.8 x 10~ 4 and 
3.7 x 10~ 4 Mq yr _1 as the convective boundary moves. The de- 
tails of the variation are described in Section |3~2l A corollary of 
the self-limiting behaviour is that the only major effect of chang- 
ing the material composition is to change its opacity. This in turn 
changes the Eddington limit and therefore the accretion rate but the 
envelope structure remains almost entirely unchanged. In the con- 
vective regions, the envelope has an adiabatic index of about 1.34 
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Table 1. Properties of the fiducial model for increasing values of Mbh- The first and last entries correspond to the initial and final models in 



the run, 


respectively. Density profiles are 


plotted in 


Fig-E 
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/10 3 K 
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0.00 


5 


2.13 


0.00 


3.48 


8.71 x 10- 


-5 


40.8 


14.3 


0.0166 


0.303 


0.51 


100 


1.79 


3.83 


2.92 


5.47 x 10" 


-8 


3.54 


5.22 


3.66 


2.09 


1.03 


200 


2.08 


24.9 


3.40 


1.30 x 10" 


-8 


2.23 


4.77 


11.1 


2.70 


2.23 


500 


2.96 


241 


4.83 


2.42 x 10" 


-9 


1.33 


4.55 


41.2 


3.54 


3.70 


1000 


3.70 


1359 


6.05 


6.49 x 10- 


10 


0.88 


4.49 


113 


4.07 


4.23 


1194 


3.53 


3360 


5.81 


3.71 x 10- 


10 


0.71 


4.51 


185 


3.96 



(corresponding to a polytropic index n « 2.90) confirming the ex- 
pectation that the envelope is approximated by an n = 3 polytrope. 
The boundaries of the convective regions depend on the ionization 
state of the gas but, for most of the evolution, all but the outermost 
few 10 Mq are convective. 

Fig. Q] shows a sequence of density profiles of the envelope 
when Mbh/Mq = 5, 100, 200, 500, 1000 and 1194 (further pa- 
rameters are listed in Table[T). These profiles demonstrate two fea- 
tures of the envelope structure. The first is the central condensation 
which can be seen from the slight rise in the density at the inner- 
most radii. This feature is clearer at smaller BH masses. It is caused 
by the lack of pressure support at the inner boundary. In order to re- 
main in hydrostatic equilibrium, the equations require 



dP 
dr 



GpM 



(16) 



at the boundary, so the pressure gradient st eepens and the density 
gradient follows. Irluntlev & Saslawl i 19751) called such structures 
loaded polytropes. 

The second feature, apparent in all but the first density pro- 
file in Fig. [T] is the density inversion in the outer layers. It ap- 
pears once the photospheric temperature T$wt drops below about 
8000 K. From then, the surface opacity increases owing to hydro- 
gen recombination. The Eddington luminosity falls and the star's 
luminosity apparently exceeds the Eddington limit. It is well known 
that hydrostatic models can susta in this super-E ddington luminos- 
ity through a density inversion dLangerlll997l) . As an additional 
check, we calculated the volume-weighted average of 37 — 4 and 
found it to be positive, indi cating dynamical (but not pulsational) 
stability jCox & Giulill 19681 Section 27.3b). 



3.2 Evolution 

The sequence of density profiles in Fig. [T] shows the interior den- 
sity decreasing over time. If we consider equation 1 1151 . ignoring 
constants and using P oc p 1 , then 



Mbh oc M, 



bhP 



(17) 



The accretion rate Mbh is held approximately constant by the Ed- 
dington limit and Mbh is always increasing. Thus, for any reason- 
able adiabatic index (7 < 3), p decreases at the inner boundary. 
Initially, the density decreases relatively rapidly and the envelope 
expands. The expansion occurs owing to the opacity peak at the sur- 
face due to hydrogen. The rate of change of the surface radius is at 
most about 0.1 Rq yr -1 which is five orders of magnitude smaller 
than the free-fall velocity. These models are thus still in hydrostatic 
equilibrium. 

Fig. [2] shows the accretion rate on to the BH as a function of 



its mass. Fig. [3] shows the locations of convective boundaries as a 
function of BH mass and demonstrates how the rapid changes of 
the BH accretion rate while Mbh < 120 Mq coincide with the 
disappearance of radiative regions. The disappearance is due to the 
decreasing density throughout the envelope. 

Before the end of the evolution the accretion rate achieves a 
local maximum. At the same time, the photospheric temperature 
reaches a local minimum and the envelope radius a maximum (see 
Fig.O. We do not have a simple explanation for this but some in- 
sight is offered in Section l4~2l Evolution beyond the results shown 
here is impossible. The code reduces the timestep below the dy- 
namical timescale indicating that it cannot construct further models 
that satisfy the structure equations. 



3.3 Termination of the run and subsequent evolution 

For the fiducial run, the evolution terminates when Mbh = 
1194Mq. The physical reason for this upper limit remains elu- 
sive but we have made some progress in understanding it using a 
modified version of the Lane-Emden equation (see Section[4). The 
existence of the limit is certainly robust as it is does not depend on 
the total mass of the quasi-star over at least two orders of magni- 
tude (see Section [5"!4t nor on whether the envelope mass changes 
in time (see Section[53}- At the end of the run the cavity contains 
a further 3360 Mq. Under our assumptions, some of this material 
is already moving towards the BH and may become part of it. If 
the BH accretes all the mass in the cavity, its final mass would 
be Mbh ~ 4554 Mq, nearly half of the total mass of the original 
quasi-star. Presuming the BH accretes at its Eddington-limited rate, 
this growth would take about 57 Myr. 

What actually happens to the material in the cavity after the 
end of the hydrostatic evolution? Because we do not model it, we 
can only speculate. It appears that the entire envelope might be 
swallowed but it is not certain that this should be the case. The ac- 
cretion flow is convective so there must be a combination of inward 
and outward flowing material within the Bondi radius. In the theo- 
retical limit for a purely con vective flow, the accretion rate is zero 
dOuataert & Gruzinovl 2000) and about half of the material is mov- 
ing inward and half outward. If the flow is sustained then we might 
expect at least half of the cavity mass to be accreted. On the other 
hand the flow structure might change completely. The infalling ma- 
terial could settle into a disc and drive disc winds or jets so that the 
overall gain in mass is relatively small. The envelope is also con- 
vective and we expect equal masses of gas to be moving inward and 
outward. Hydrostatic equilibrium is presumably failin g so the dy- 
namic s might change drastically here, too. Recently. I Johnson et al.l 
(201 1) modeled the accretion on to massive BHs formed through 
direct collapse. They assumed that the BH accretes from a multi- 
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Figure 2. Plot of BH accretion rate A/bh against BH mass AfgH for various 
pairs of values of the radiative and convective efficiencies (e, r]). The fidu- 
cial values (0.1, 0.1) correspond to the solid line. As expected, changing 
the radiative efficiency changes ATbh but leaves the overall structure unaf- 
fected. Decreasing jj causes the envelope to be hotter and denser in order to 
achieve the same luminosity shifting the discontinuities in the gradient of 
the accretion rate to later times. 
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Figure 3. Radial location of convective boundaries and the Bondi radius 
for the fiducial run. The shaded regions are the radiative parts of the en- 
velope. The outer solid lines show the extent of the hydrostatic envelope. 
The qualitative change in the locations of convective boundaries causes the 
discontinuities in the gradient of the accretion rate that are seen in Fig. [2] 
The outermost layer of the envelope is radiative (as it should be) but it is 
too narrow to be seen here. 



colour black body disc after its quasi-star phase and found that, 
once the BH mass exceeds about 10 4 Mq, the accretion rate de- 
creases due to radiative feedback. This result supports the case for 
a substantial decrease in the BH growth rate if a thin disc forms af- 
ter the quasi-star phase but additional growth can occur during the 
transition to a new structure. 



4 FINAL BH MASS LIMIT 

We have not yet found a simple explanation for the existence of 
the upper limit to the BH mass in a quasi-star. In this section, we 
explore some theory that confirms the existence of such a limit and 
points towards possible mechanisms. 



4.1 Loaded polytropes 

Consider the equations of hydrostatic equilibrium and mass conser- 
vation truncated at some radius ro and loaded with some mass Mo 
interior to that point. The equations are 



dP _ Gmp(r) 
dr r 2 
and 

= Anr 2 p(r) 

dr 



(18) 



(19) 



with the central boundary conditions m\ ro = Mo, where ro and 
Mo are fixed. We scale the pressure and density using the usual 
polytropic assumptions 



P = Kp 1+ ~ 
and 

p = po9". 

We define the dimensionless radius by 

r = af, 

where 

2 (n+l)K i-i 
° = —i^ pS ■ 

We scale the mass interior to a sphere of radius r by defining 

m 



(20) 



(21) 



(22) 



(23) 



(24) 



</>(£)= . 

47rpoa 

dHunflev & Saslawll 19751) . The non-dimensional form of the equa 
tions is then 



and 



(25) 



(26) 



with boundary conditions 0(£o) = 1 and 0(£o) = 4>o (where, by 
definition, £o = m/oOE] 



1 If one takes go = <t>0 = 0, differentiates equation (25j and substitutes for 
using equation (26) , one arrives at the usual Lane-Emden equation. 
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Figure 4. Plot of the ratio of BH mass against total mass against inner enve- 
lope radius to outer envelope radius. The double-valued curves correspond 
to polytropic solutions (with n = 2, 3 and 4). The monotonic curve (solid 
line) represents the fiducial results. The polytropic models show the exis- 
tence of an upper limit to the BH mass ratio. The fiducial results reach a 
similar limit, but the mass ratio cannot decrease, so the evolution termi- 
nates. 



We can express 0bh (the scaled BH mass) in terms of the inner 
radius by rescaling equation (|2j as follows. 



r 2G Mbh 
so = — = 5— 

a a c% 

2G a 3, n 
= in poa 0BH — 

(n + l)Kp S 

o 2, 4ttG 
= 2na </>bh ; — - 



(n + l)KpS 



2n^BH- 



(27) 
(28) 

(29) 
(30) 



Similarly, we can derive the following relation for Af cav from equa- 
tion {3). 



47rpoa 3 



- V 

-g€o. 



(3D 
(32) 



The scaled mass and radius boundary conditions are now related by 



0o = <K£o) 



1 2 3 

—to + -So- 

2n 3 



(33) 



Thus, for a given polytropic index n, we can choose a value £0 an d 
integrate the equations. 

We integrated a sequence of solutions for n = 3 and found 
that there exists a maximum value of 4>bh/4 i * — 0.102, where 0* is 
the scaled total mass of the quasi-star. The maximum occurs when 
£0 = 0.995. We investigated polytropic indices between 2 and 4 
and found a similar limit in all cases. For n = 2 we found a max- 
imum </>bh/<^* = 0.127 when £0 = 1.012 and for n = 4 the 
maximum was 4>bh/4>* = 0.089 at £0 = 0.968. Fig.|4]shows plots 



of the curves of the ratio of inner to outer envelope radius (£0 /£* 
in scaled variables, where is the outer radius of the envelope) 
against the BH fractional mass (0bh/0* in scaled variables) for 
n — 2,3 and 4 together with our results for the fiducial model. The 
maximum mass ratio is clear in each curve. In principle, further hy- 
drostatic solutions probably exist along the sequence computed by 
STARS but they require that the BH mass decreases. 

The conclusion we draw is that the models produced by the 
code come to a halt because no further realistic hydrostatic solu- 
tions can be found along the sequence. The discussion in this sec- 
tion, even though it does not immediately offer an explanation of 
why, thus indicates that such a limit does truly exist. The limit is 
robust in the sense that it exists at approximately the same ratio for 
all masses and does not depend on details such as opacity or energy 
transport (as long as the envelope is convective). The main depen- 
dence in real models is due to variation in the polytropic index. 



4.2 The U-V plane 

It is known that (see Chandrasekhail ll939L Section 4.8), if one 
has determined a solution (?(£) to the Lane-Emden equation, then 
#'(£') = C ™ _1 0(C£) is also a solution for an arbitrary constant 
C. By selecting appropriate variables that are invariant to such a 
transformation, one can find all the related solutions in a single 
calculation. Tho ugh a number of appropriate variables exist, many 
authors (such as lKimura|[l98l1 ; iHoredJl 1 9871) use U and V, defined 
below. Converting to the new variables also allows the analysis of 
solutions of the polytropic equation that do not extend to £ = 
or contain the entire mass of the envelope. Physical solutions are 
restricted to positive U and V. The first quadrant of the plane that 
they form is a useful tool for analyzing the behaviour of solutions. 
We begin by defining 



U = 
and 

V = 



(34) 



(35) 



Let us differentiate the logarithms of these variables with respect to 



3 


n d9 


1 d(f> 


! + 


d£ " 


~ 0< 


3 


n(j> 




t ~ 


6^ ~ 




5< 3 


-nV 


-U) 



and 

1 dV 
V~dJ 



ld9 \d(j> 

0df + ~d£ 



1 ±_ eo n 

s" <t> 

1{-1 + U + V). 



(36) 
(37) 
(38) 

(39) 
(40) 
(41) 



Dividing these two equations eliminates £ and allows us to write 
the first-order equation 



dV_ 

dU 



v_ f u + v-i 

U \ U + nV - 3 



(42) 
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Figure 5. Some features of the U—V plane forn = 3. The dashed lines cor- 
respond to limiting values of dV/dlf, as indicated on the plot. The dotted 
line shows the curve of central boundary conditions for the loaded poly- 
tropes described in Section |4~71 with BH mass increasing towards the top 
and right. The solid line shows the standard n = 3 polytrope. The dash- 
dotted line is the loaded polytrope corresponding to the largest BH-envelope 
mass ratio. The intersection of the dU / dV = line and the curve of initial 
values corresponds well with the maximum accretion rate in in Fig. [2] It is 
interesting that the final model takes initial values approximately halfway 
between this point and the intersection with the n = 3 polytrope. 



Polytropes that have zero mass and non-zero density at £ = begin 
at(tr,V) = (3,0). 

Using the scaled boundary conditions in the previous section, 
we can determine the curve in the U—V plane which corresponds 
to the interior values of our models. The relevant initial points are 
given by 



U = 



-e 2 1 + V 



and 



Vo 



9o 



1 2 2 
2^ + 3 €o ' 



(43) 



(44) 



We can eliminate £o to give the curve of initial points in the U—V 
plane. 



Vo = — (3-2(7o) _1 . 
in 



(45) 



Different BH masses correspond to different points along this 
curve. A solution for a particular BH mass is then a curve going 
from the relevant point on the (Co, Vb) curve towards (7 = 0. 

Fig.|5]shows a number of features in the U -V plane for n = 3. 
Along the dV/dU — or dU/dV = lines, the solutions are hor- 
izontal or vertical, respectively, in the U-V plane. All solutions 
bounded below by the dU jdV = have decreasing U and increas- 
ing V everywhere. 

How do the points of interest correspond to our models? The 
intersection of the dU/dV — line and the (Uo, Vo) curve cor- 
responds to £n = 0.835, for which the BH-envelope mass ratio is 



0.0953. This appears to correspond to the point at which the ac- 
cretion rate achieves its final local maximum. The dash-dotted line 
indicates the solution to the equations in Section |4~T| that has the 
largest BH-envelope mass ratio. The curve's initial point appears to 
be about halfway between the intersections of the (Uq,Vo) curve 
with the dU I dV = and with the polytropic solution. Though the 
mechanism behind the mass limit is not yet known, the analysis 
in this section shows that the limit is reproduced in the polytropic 
approximation, so the relevant physics is contained within mass 
conservation and hydrostatic equilibrium. 



5 FURTHER RESULTS 

In Section[3] we established the basic qualitative structure and evo- 
lution of the quasi-star envelope. In this section, we explore their 
dependence on some of the parameters of the model. 



5.1 Radiative and convective efficiencies 

We first vary the accretion rate by adjusting the parameters e and 
r\. Fig.[2]shows the accretion history against BH mass for a number 
of parameter choices. Because the luminosity always settles on the 
same convection-limited rate, changing e only rescales the accre- 
tion rate through equation ((9) and has no effect on the structure. 
The final BH mass and intermediate properties are the same. The 
only difference is that the evolution takes longer for larger values 
of e. 

Increasing the convective efficiency r\ allows a greater flux to 
be radiated. Because it is a limiting factor, a larger value of r\ al- 
lows a larger accretion rate for given interior conditions. In order 
to establish the same overall luminosity, the envelope must be less 
dense. This explains the different times at which the discontinuities 
in the gradient of the accretion rate appear. 



5.2 Cavity properties 

The second property we adjust is the radial profile of the density 
inside the cavity. If the inward transport of angular momentum by 
advection and convection is greater than the outward transport by 
magnetic fields and other sources of viscosity, then the density pro- 
file tends towards p(f) oc r~ 5 . Recomputing the cavity mass from 
equation ([3} gives 



M cav = —p r B . 
5 



(46) 



The quasi-star evolution is shown in Fig. [6] The smaller interior 
mass leads to a less evident density spike at all times. The decreased 
interior mass is subject to a lower mass limit for the BH. The final 
BH mass is 1077 Mq. 

Although the change to the cavity properties affects the nu- 
merical results, there is no qualitative change to the quasi-star evo- 
lution. The evolution terminates owing to the same physical mech- 
anism as is analysed in Section|4] 



5.3 Envelope mass loss and gain 

To illustrate the effect of a net accretion rate on to the surface of 
the envelope, Fig. [6] shows the evolution of the fiducial run if the 
envelope is accreting at a constant rate of 2 x 10~ 3 Mq yr _1 . Al- 
though this rate initially exceeds the quasi-star's Eddington limit, 
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Figure 6. Plot of BH accretion rate A/bh against BH mass Mbh for the 
fiducial ran, a ran w ith constant accr etion on to the surface of the star ("Env 
Acc"), a ran with a Reimers ( 1975) mass loss rate ("Reimers") and a ran 
with a shallower radial dependence of the interior density ("p oc r~ V 2 "). 
The Reimers rate leaves an envelope of about 7750 Mq and the BH is 
proportionally smaller. The accreting envelope also leaves a proportionally 
scaled BH. A shallower radial dependence of the interior density leads to a 
smaller cavity mass and therefore a smaller interior density and accretion 
rate. 



such rapid infall is believed to occur as long as the bars-within- 
bars mechanism is transporting material towards the centre of the 
pregalactic cloud. Once the surface temperature of the quasi-star 
decreases below about 8000 K, the Eddington limit becomes much 
larger owing to the decreasing opacity. The accreted mass is simply 
added to the surface value of the mass co-ordinate and no additions 
are made to any other equations. In particular, we do not include 
a ram pressure at the surface. The only qualitative change to the 
evolution is that it takes longer than if the total mass had been held 
constant at the same final value of 17 540 Mq. The final BH mass 
is subject to the same ratio limit so a larger final envelope permits 
a larger final BH. 

To investi gate the effect of mass loss we use a Reimers rate 
jReimersll 1975b . It is an empirical relation to describe mass loss in 
red giants. The mass-loss rate is 



M, oss = 4 x 10" 



;L*R, M(T 



© 



Ah 



(47) 



L R - 

Fig. [6] shows the evolution for quasi-star envelopes with this pre- 
scription. The mass loss is significant but, again, no qualitative 
change in the results is seen. The limit holds and the BH is slightly 
smaller, precisely in proportion with the decrease in the mass of the 
envelope. 



5.4 Initial envelope mass 

The structure of the envelope appears to be chiefly dependent on the 
ratio of envelope mass to BH mass. We constructed a set of models 
with initial total masses A'h / Mq = 10 4 , 2 x 10 4 , 5 x 10 4 , 10 5 , 2 x 



Figure 7. Plot of the evolution of quasi-stars of different total masses. The 
BH mass and accretion rates have been divided by the total quasi-star mass 
to illustrate the consistency of the upper mass ratio limit of 0.119 and the 
slight dependence of the accretion rate with quasi-star mass. Because larger 
quasi-stars permit greater scaled accretion rates, they have shorter hydro- 
static lifetimes. 



10 5 , 5 x 10 5 , and 10 6 to confirm this and explore how the temporal 
evolution of the envelopes is affected by total mass. 

Fig. [7] shows the evolution of the various envelopes with the 
BH mass and accretion rate divided by the quasi-star masses. We 
note that the fractional upper BH mass limit holds for all the quasi- 
star masses in this range. Also, the temperature profiles by mass of 
the envelopes depend almost exclusively on the BH-envelope mass 
ratio. The variation in inner or surface temperature, for a given mass 
ratio, with respect to the total mass of the quasi-star is less than 0.05 
per factor of ten in the total mass. 

The density and radius profiles are more strongly dependent 
on the mass of the envelope. We find that, for a given mass ratio 
and once the entire envelope has become convective, the following 
approximate relations hold for the properties of two quasi-stars of 
different masses (denoted by subscripts 1 and 2). 



R.,i 

R*,2 

po,i 

P0,2 

and 



Ah,i 
~Ah~~2 



M,,i 
~Ah~2 



BH.l 



\MbH,2 



M,,i 



(48) 



(49) 



(50) 



For example, compared to a quasi-star of mass 10 4 M at the same 
BH-envelope mass ratio, a 10 5 Mq quasi-star will have an outer 
radius that is 10 ' 54 = 3.47 times greater, an interior density that 
is lO 66 = 4.57 times smaller and a mass accretion rate that is 
10 114 = 13.8 times greater. 

The final relation implies that the lifetime of a quasi-star scales 
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Figure 8. Comparison of analytic estimates of BRA08 (dashed lines) 
against results for our fiducial run (solid lines) for interior temperature (top), 
surface temperature (middle) and envelope radius (bottom). BRA08's esti- 
mate of the interior temperature is accurate but those for the photospheric 
temperature and envelope radius become increasingly inaccurate as the BH 
grows. 



as tqs oc so larger quasi-stars have slightly shorter hydro- 

static lifetimes. Fig. [7] shows how the BH-envelope mass ratio for 
which the entire envelope is convective also depends on the enve- 
lope mass. This has a small effect on the lifetime of the quasi-stars. 
By fitting a straight line to the log TQs-log M* relation for the seven 
models here, we find that the lifetimes scale as tqs oc M^ 013 . 
More precisely, we find 

Iogio tqs = -0.126 log 10 M, + 7.12 (51) 

Note that M» here denotes the initial mass of the quasi-star. In all 
other relations M* slowly decreases during the quasi-star's evolu- 
tion owing to the mass-energy lost as radiation. 



5.5 Comparison with Begelman et al. (2008) 

BRA08 have estimated some envelope properties presuming that 
the envelope is described by an n = 3 polytrope. Their estimates 
are made in terms of an overall accretion efficiency paramete^] 
Qbra which is determined by numerical factors in the accretion rate 
including the radiative efficiency, convective efficiency and adia- 
batic index. Because our adiabatic index is not fixed, we have cal- 
ibrated qbra using the BH luminosity (equation 3 of BRA08). To 
compare the fiducial run, we selected qbra ~ 0.257 which makes 
their analytical BH luminosity accurate to within 0.2 per cent over 
the entire evolution. 

We compare the following estimates for the inner temperature, 
photospheric temperature and envelope radius (equations 7, 11 and 



2 This should not be confused with a defined in Section l4~Tl 
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Figure 9. Plot of BH accretion rate \l\w\ against BH mass Mbh for dif- 
ferent choices of the inner boundary radius. The various values of ro lead 
to qualitatively similar results but the quantitative evolution is strongly af- 
fected. 

10 of BRA08, respectively). 

Tai = 1.0 X 10 — ctBRA-r-rr- T! - K 

(53) 

and 

(^)*~ 

(54) 

Here, Lem = 4-kGcM/k is the Eddington luminosity. BRA08 
computed this using the opacity at the boundary of the convective 
zone but such estimates differ by a factor of the order of k/k cs 
when compared with our results. Our comparison is made using 
the Eddington limit with opacity Kes = 0.34 cm 2 g _1 . 

In Fig. [8] we plot the three estimates against the results from 
our fiducial run. The estimate for the interior temperature is accu- 
rate to within 20 per cent. The deviation grows as the approxima- 
tion of the envelope to an n = 3 polytrope becomes increasingly 
less accurate. 

The estimate for the photospheric temperature is highly inac- 
curate. At the end of the run the estimated photospheric temperature 
is about 2400 K compared to the model result of about 4500 K. Be- 
cause the BH luminosity estimate is very accurate, it then follows 
that the envelope radius is inaccurate. The surface luminosity must 
be related by L« = 7raci? 2 Tj f . This is confirmed in the bottom 
panel of Fig. [8] 

BRA08 argue that quasi-star evolution terminates owing to the 
opacity at the edge of the convection zone increasing. The increased 
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Figure 10. Contours of (j>o/<p* for the models calculated for each pair of 
initial conditions (£o , <t>0 ) in the plane. There is a clear boundary along 
<t>Q = £o near the origin. Solutions above this line must have negligibly 
small envelopes by mass. 



Figure 1 1. Contours of 0o I4>* for the models calculated for each pair of ini- 
tial conditions (Uo, Vo) in the plane. The boundary in Fi g . 1 1 I corresponds 
to the critical point (0,1). 



opacity causes the envelope to expand and the opacity increases fur- 
ther. The envelope then expands further and the process is claimed 
to run away. BRA08 refer to this process as the opacity crisis. Our 
results do not terminate for this reason. Similar behaviour does oc- 
cur at the beginning of the evolution while the photospheric tem- 
perature is greater than 10 4 K but it does not disperse the quasi-star. 
For most of a quasi-star's evolution, the opacity at the convective 
boundary is already beyond the H-ionization peak and is decreasing 
as the BH grows. 



6 SENSITIVITY TO THE INNER BOUNDARY 

Besides the physical parameters of the models, we have also con- 
sidered the effect of changing the inner boundary radius. Fig. [9] 
shows the evolution of quasi-stars where the inner radius was 
changed to half and twice the Bondi radius. The final BH mass 
is strongly affected although the evolution is qualitatively similar. 
The results of our polytropic analysis in Section|4]are also affected 
in a consistent manner. Although reasonable, our choice of inner ra- 
dius is somewhat arbitrary and critical in deciding the quantitative 
evolution of the quasi-star. A self-consistent model to determine 
the appropriate value of ro is thus highly desirable but none is yet 
available. We continue to seek a suitable resolution to this issue. 

The sensitivity of the models to the inner boundary precipi- 
tated further interesting results. We note that the final BH mass is 
scaled by approximately the same factor as the inner radius. That 
is, if ro is doubled, the BH reaches about half the final mass that 
it reached before. If the material properties at the inner boundary 
were largely unchanged, this would imply that the inner radius at 
the termination of the runs is approximately the same. The output 
from the models indicates that ro / Rq in the final reliable model 
of each run in Fig.|9]is approximately 2.05 x 10 4 , 1.94 x 10 4 and 
1.79 x 10 4 for ro = 2rs, rs and rB/2. This indicates that the limit 



determined for the polytropic models in Section |4]may be related 
more strongly to the ratio of inner to outer radius than the ratio of 
interior to total mass. 

It may appear from Fig. [4] that the limiting radius ratios for 
the STARS output and the polytropic models are inconsistent. The 
discrepancy is due to the ionization zones inside the envelope. It 
is also difficult to identify the appropriate values because the in- 
ner and outer radii are changing most rapidly at the end of the 
evolution. Using the last models in each run of Section [5T41 that 
have converged in thermal equilibrium, the limit is approximately 
r /R* w 0.460. 

We additionally found that we could not construct model en- 
velopes for ro < 0.3rB and find that this is reflected in the poly- 
tropic models. To construct Fig. [10] we have calculated the ratio 
4>o /(f), for each (£o> 4>o) pair in the plane and plotted contours from 
to 1 in steps of 0.1. Near the origin, there is a clear boundary 
along the line 4>o — Co- The boundary softens for larger values of 
4>o and Co- For small Co, models above the boundary necessarily 
have 0o/0* ~ 1. In other words, any model with a sufficiently 
small inner boundary must have an appropriately small envelope. 

The limit becomes clearer when plotted in the U -V plane dis- 
cussed in Sectionf4j2] We have done so in Fig.QT] The contours are 
divided across the point (0,1). This is a critical point in the U—V 
plane for all polytropic indices. It is neutrally stable in the direction 
(1,-1) and unstable along the V-axis i.e. (0,1). 

The present analysis thus provides several important results. 
First, the choice of inner radius is crucial to determining the evolu- 
tion of the envelopes. The Bondi radius is a reasonable choice but 
there is no evidence to confirm that it is the correct one. Secondly, 
the limit determined in Section [4] seems more strongly related to 
the radius ratio than the mass ratio. This was not obvious at first 
because there was no apparent reason to modify the inner radius. 
Thirdly, there is a limit to the smallest value that can be chosen for 
the inner radius. 
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7 CONCLUSION 

We have modified the Cambridge STARS stellar evolution code to 
model the evolution of quasi-star envelopes. Our first new result is 
the existence of a robust upper limit on the ratio of inner BH mass 
to the total mass, equal to about 0.119, of the system in hydrostatic 
equilibrium. The limit is reflected in solutions of the Lane-Emden 
equation, modified for the presence of a point mass interior to some 
specific boundaries. After considering variation of the inner radius, 
this limit is possibly better interpreted as a limit to the ratio of inner 
radius to envelope radius. The value of the limit from the STARS 
output is about 0.46. 

All the evolutionary runs here terminate once the limit is 
reached. It is difficult to say what happens to the BH and envelope 
after the hydrostatic evolution ends. Some of the material within the 
Bondi radius has begun accelerating towards the BH so we expect 
that it can be captured by the BH. The remaining material may be 
accreted or expelled, depending on the liberation of energy from the 
material that does fall inwards. After the BH has evolved through 
the quasi-star phase, it is probably limited to accreting at less than 
the Eddington limit for the BH. 

The models presented are crucially sensitive to the choice of 
inner boundary radius and the results should be treated with due 
caution. While the Bondi radius used here is reasonable, we con- 
tinue to seek a less arbitrary set of boundary conditions. 

In light of these results, it appears that quasi-stars produce 
BHs that are on the order of at least 0.1 of the mass of the quasi-star 
and around 0.5 if all the material within the inner radius is accreted. 
For conservative parameters, this growth occurs within a few mil- 
lion years after the BH initially forms. Realistic variations in the 
parameters (e.g. larger initial mass, lower radiative efficiency) lead 
to shorter lifetimes. Such BHs could easily reach masses exceed- 
ing 10 9 M early enough in the Universe to power high-redshift 
quasars. 



ACKNOWLEDGEMENTS 

The authors would like to thank Peter Goldreich for valuable 
comments. We are grateful to Mitch Begelman, particularly for 
comments that led to the investigation in Section [6] CAT thanks 
Churchill College for a Fellowship. 



REFERENCES 

Appenzeller I., Fricke K., 1971, A&A, 12, 488 
Begelman M. C, 1978, MNRAS, 184, 53 
Begelman M. C, 2010, MNRAS, 402, 673 
Begelman M. C, Shlosman I., 2009, ApJL, 702, L5 
Begelman M. C, Volonteri M, Rees M. J., 2006, MNRAS, 370, 289 
Begelman M. C, Rossi E. M., Armitage P. J., 2008, MNRAS, 387, 1649 
(BRA08) 

Bond J. R., Arnett W. D., Can- B. J., 1984, ApJ, 280, 825 

Bondi H., 1952, MNRAS, 112, 195 

Bromm V., Larson R. B., 2004, ARA&A, 42, 79 



Chandrasekhar S., 1939, An Introduction to the Study of Stellar Structure. 
Chicago Univ. Press, Chicago 

Clark P. C, Glover S. C. O., Klessen R. S., 2008, ApJ, 672, 757 

Cox J. P., Giuli R. T., 1968, Principles of Stellar Structure. Gor- 
don & Breach, New York 

Devecchi B., Volonteri M., 2009, ApJ, 694, 302 

Dijkstra M., Haiman Z., Mesinger A., Wyithe J. S. B., 2008, MNRAS, 

391, 1961 
Eddington A. S., 1918, ApJ, 48, 205 
Eggleton P. P., 1971, MNRAS, 151, 351 
Eggleton P. P., 1972, MNRAS, 156, 361 
Eggleton P. P., 1973, MNRAS, 163, 279 

Eggleton P. P., Faulkner J., Flannery B. P., 1973, A&A, 23, 325 
Eldridge J. J., Tout C. A., 2004, MNRAS, 348, 201 
Fan X. et al., 2006, AJ, 131, 1203 
Flammang R. A., 1984, MNRAS, 206, 589 
Fowler W. A., 1964, Rev. Modern Phys., 36, 545 

Freese K., Bodenheimer P., Spolyar D., Gondolo P., 2008, ApJL, 685, 
L101 

Fryer C. L., Woosley S. E., Heger A., 2001, ApJ, 550, 372 
Horedt G. P., 1987, A&A, 177, 117 
Huntley J. M., Saslaw W. C, 1975, ApJ, 199, 328 
Itoh N, Hayashi H., Nishikawa A., Kohyama Y., 1996, ApJS, 102, 411 
Jiang L. et al., 2008, AJ, 135. 1057 
Johnson J. L., Bromm V., 2007, MNRAS, 374, 1557 
Johnson J. L., Khochfar S., Greif T. H., Durier F, 2011, MNRAS, 410, 
919 

Kimura H., 1981, PASJ, 33, 273 

King A. R., Pringle J. E., 2006, MNRAS, 373, L90 

Langer N, 1997, in Nota A., Lamers H., eds, ASP Conf. Ser., Vol. 120, 

Luminous Blue Variables: Massive Stars in Transition, p. 83 
Lodato G., Natarajan P., 2006, MNRAS, 371, 1813 
Loeb A., Rasio F. A., 1994, ApJ, 432, 52 

Milosavljevic M., Couch S. M., Bromm V., 2009, ApJL, 696, L146 
Narayan R.. Igumenshchev I. V, Abramowicz M. A., 2000, ApJ, 539, 798 
Ohkubo T., Umeda H., Maeda K., Nomoto K., Suzuki T., Tsuruta S., Rees 

M. J., 2006, ApJ, 645, 1352 
Pols O. R., Tout C. A., Eggleton P. P., Han Z., 1995, MNRAS, 274, 964 
Quataert E., Gruzinov A., 2000, ApJ, 539, 809 
Regan J. A., Haehnelt M. G., 2009, MNRAS, 393, 858 
Reimers D., 1975, Mem. Soc. R. Sci. Liege, 8, 369 
Schleicher D. R. G., Spaans M., Glover S. C. O., 2010, ApJL, 712, L69 
Shang C, Bryan G. L., Haiman Z., 2010, MNRAS, 402, 1249 
Shapiro S. L., 2004, in Ho, L. C, ed., Carnegie Observatories Astroph. 

Ser., Coevolution of Black Holes and Galaxies. Cambridge Univ. Press, 

Cambridge, p. 103 
Shlosman I., Frank J., Begelman M. C, 1989, Nat, 338, 45 
Spaans M., Silk J., 2006, ApJ, 652, 902 

Spolyar D., Freese K., Gondolo P., 2008, Phys. Rev. Lett., 100, 051101 
Tegmark M., Silk J., Rees M. J., Blanchard A., Abel T., Palla F, 1997, 
ApJ, 474, 1 

Thorne K. S., Zytkow A. N, 1977, ApJ, 212, 832 

Volonteri M., 2010, A&AR, 18, 279 

Willott C. J. et al., 2010, AJ, 139, 906 

Wise J. H, Turk M. J., Abel T., 2008, ApJ, 682, 745 

This paper has been typeset from a TpX/ KTpX file prepared by the 
author. 



© 2011 RAS, MNRAS 414, 2751-2762 



